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Abstract. The matrix equation XA + AX T = 0, which has relevance to the study 
of Lie algebras, was recently studied by De Teran and Dopico (Linear Algebra 
Appl. 434 (2011), 44—67). They reduced the study of this equation to several 
special cases and produced explicit solutions in most instances. In this note we 
obtain an explicit solution in one of the difficult cases, for which only the dimen- 
sion of the solution space and an algorithm to find a basis of this space were 
known previously. 



1. Introduction 

In the recent article [1J, F. De Teran and EM. Dopico considered the equation 

XA + AX T = 0, (1) 

where A is a given square complex matrix. For each fixed A, the solution space 
to CD is a Lie algebra q(A) equipped with Lie bracket [X,Y] := XY — YX. To 
characterize, up to similarity, all possible Lie algebras g(A) which can arise in this 
manner, it turns out that one need only consider those A which are in congruence 
canonical form (recall that two matrices A and B are congruent if A = P T BP for 
some nonsingular P) [TJ Sect. 2]. 

To be more specific, it is known that each square complex matrix is congruent 
to a direct sum, uniquely determined up to a permutation of summands, of three 
types of canonical matrices (see |5, Thm. 2.1] for complete details). In general, the 
Lie algebras corresponding to these canonical matrices are far from semisimple 
and many of them display quite striking visual patterns (e.g., [1 , p. 52, p. 54]). 

We are concerned here with canonical matrices of Type II. These are the matri- 
ces 

HM:= \jM o. ' 

where l n denotes the n x n identity matrix and }„(}i) is the n x n Jordan matrix 
whose eigenvalue u satisfies ^ u ^ ( — 1)" +1 . De Teran and Dopico were unable 
to find an explicit description of the Lie algebra fl(H2„(( — 1)")), although they 
computed the dimension of this algebra for even n using an algorithmic approach 
and stated the corresponding results for odd n without proof [1. Appendix A]. 
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We aim here to provide an explicit solution to the matrix equation 

XH 2 „( F ) + H 2n (ii)X T = 0. (2) 

In particular, our method does not significantly depend upon the parity of n. 
Besides providing an explicit description of the Lie algebra g(H2«((— 1)")), our 
approach is also notable for its use of a family of matrices derived from Pascal's 
Triangle. 



2. Solution 

We begin our consideration of the matrix equation <(2j by partitioning X con- 
formally with the decomposition of H2„(}i): 

'A B 
C D ' 



X 



Substituting the preceding into 10 yields the system of equations 

B} n ( F )+B T = 0, (3) 

A + D T = 0, (4) 

DJ n (ii) + J n ( } i)A T = 0, (5) 

C + J n (ji)C T = 0, (6) 

which together are equivalent to 101. Before proceeding, let us note that solving 
for X in (0, taking the transpose of the resulting expression, and substituting the 
final result back into 10 shows that X commutes with the matrix 

H 2 „( F )H 2 - T ( F ) = 7„0<) © J- T ( F ). (7) 

As noted in (TJ Lem. 6], these observations are sufficient to handle the case y. 7^ 
±1. In this case, the matrices J^ T {}i) and J n (}i) have distinct eigenvalues whence 
B = C = by [4, Cor. 9.1.2]. Since D commutes with } n (pi) by 10 and 10, it 
follows from H Thm. 9.1.1] that D is an upper-triangular Toeplitz matrix and 
that X = ( — D T ) © D. 

Now suppose that \i = ( — 1)". As before, we see that D is an upper- triangular 
Toeplitz matrix and that A = —D T . The entire difference between our two cases 
lies in the structure of the off-diagonal blocks B and C. 

We first focus our attention upon the matrix B. Taking the transpose of 10 and 
substituting the result back into 10 yields 

B/«00 = In T (F)B- (8) 
Since ]i = ( — 1)", it turns out that the Jordan canonical form of }~ T (}i) is precisely 
Jn(}i)- In fact, the similarity is implemented by a special matrix derived from 
Pascal's Triangle. 



Lemma 1. We have 



(9) 



ON THE MATRIX EQUATION XA + AX T = 



3 



where A„ is the n x n matrix whose entries are given by 



n — i 



(10) 



Zn particular, we observe the convention that (?) = whenever k > n. Thus for n 
1, 3, 5, . . . roe ziaoe 



Ai 



-1' 
12 
-1 -1 -1 



[1]/ A 3 = 
and so forth. For n = 2,4,6, .. . we have 



A 5 



-1" 

1 4 

0-1-3 -6 

12 3 4 

-1 -1 -1 -1 -1 



A 2 



-1 

1 -1 



, A 4 



-1' 
1-3 

0- 12-3 

1- 11-1 



A fi = 



"0 














-1" 














1 


-5 











-1 


4 


-10 








1 


-3 


6 


-10 





-1 


2 


-3 


4 


-5 


1 


-1 


1 


-1 


1 


-1 



and so on. 



Pf.of Lemma\T\ We first note that 10 is equivalent to 

An = /„ T (F)A„/n(f<). (11) 

Upon writing 

/«(?') = F^n + 7»(0)/ 

we find that fTTjl is equivalent to 

/„ T (0)A„/„(0) + hA„/„(0) + F/J(0)A fi = 0. (12) 

Using Pascal's Rule and the fact that }i — ( — 1)", the (i,j) entry in the right hand 
side of l fl2l l is 

[A„] + u[A„] i/hl + u[A n ]i-l,j 



0. 



7-2 
n — i + 

7-2 
n - z + 1 



7-1 



7-2 
n — i 



7-1 
n-i + 1 



(13) 



For t + j < n, we have used the fact that all of the binomial coefficients appearing 
in (fl3b vanish. □ 

Lemma 2. The inverse of A n is the reflection of A n with respect to its center. In other 
words, 

[An" 1 ]/,; = [A] w _ f+1/B _y +1 

for 1 < i,j < n. 

Pf. ofLemma^ Let A„ denote the reflection of A„ with respect to its center. First 
note that [A n A n ]y = holds whenever i > j, due to the "triangular" structure of 
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the matrices involved. On the other hand, if j > i, then it follows that 

n 

[A„A„], V = Y^[An]i,k[&]k,j 
k=\ 
n 

= Y^\.^n-i+l,n-k+l[A]k,j 
k=l 



' i+ ' ( I 1 (- i )'- <M, G- 1 )(^ 1 



t=i-\ V 1 



i + 1 



/■ i\ /— i /■ 



^ +/ (- 1 )' +/ ({_l)E(- 1 ) r 



r 



Thus A„A„ = I„, as claimed. 



Lemma 3. TTze identity 

(A„/^( F )) T = -A n /r^( H ) 

holds for £ — 0,1, . . . ,n. 



Pf. ofLemma\3\ We first claim that it suffices to prove that 

holds. Indeed, using © and (TT5l it follows that 

A£ = -A w /£(p) 

= -(A„/^( F )A- 1 )A„/r f ( F ) 

= -(/2oor r MT'oo, 

which implies <fl4]> . 

In order to prove 1151 1, it suffices to derive the equivalent identity 
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Computing the entry in the matrix product A~ 1 Aj, we see that 

k=l 
n 

= Yli[ A n\n-i+\,n-k+l[^n]j,k 
k=l 



j2 ly{-Y-?~ l y n _j 



k=i 

'n - k\ (k - 1 



i+j+l 



1=0 



i — 1 J \n — j 



-"• +i+ \n + "i-i) < 17 » 



where the passage from fhSt to lUZjl follows by computing the coefficient of x n_1 
in the identity 

x n-j 



(l-x)"-/ +1 (l-x) n +H+i' 



Having digressed somewhat upon the properties of the matrix A„, let us return 
now to the equation ^8). Upon substituting © into ®, we find that 

(A" 1 B)/„( ? /)=/„( F )(A„- 1 B) / 

from which it follows that B = A n U for some upper triangular Toeplitz matrix 
U 2J Thm. 9.1.1]. However, it is important to note that not every such U leads to 
a solution of the original equation (O. 

Now observe that the matrices I n , J n (pi), . . . , (ji) form a basis for the space 
of n x n upper triangular Toeplitz matrices. To see this, simply note that J*(ji) 
contains nonzero entries on the kth superdiagonal whereas the fcth superdiago- 
nals of the matrices l n , J n (ji), • • •/ /h _1 (f) are identically zero. Since U is an 

upper triangular Toeplitz matrix, we may write 

;=o 

for some constants ocq, ol\, . . . , ot n _\. 
Substituting B — A„U into © yields 

A„U=-(A„U) T ]- 1 ( } i) r (19) 
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which, in light of dl8t , provides us with 



n-l 



n-l 



£ Kj&nj'n(F) = - E a /c(An/£(F)) In V)- 
;=0 k=0 

Using the identity (TT4T ) in the preceding formula, we obtain 



;=0 



jt=0 



which can be rewritten as 



n-l 



E(«f -«w-l-*)/nO<) = °- 
1=0 

Since the matrices I n , J n (}i), J%(}i), . . are linearly independent, it follows 

that the matrix i flBI l is determined by exactly [£] free parameters. Thus the gen- 
eral solution to © is 




For « = 3, the general form of the matrix B is 

2(«i -«o)l [ -ft' 

2(«o-«i) 2(a -«i) = jS a ft 
_2(«i-«o) -«o J [-ft — jS . 

where /So = «o an d ft = 2(ao — For n = A we have 







_2a + 2«j 
and for n = 5 we obtain 



-2a - 2a! 

2a + 2«i — 3«o — 3aj 

-2«o — 2«i HQ + «i — 3ao — oc\ 

kq + oti 2a — «o 



o 
o 


2a - 2«i + 2a 2 



-2a + 2ai - 2a 2 2a - 2«j + 2« 2 



-2a + 2a! - 2a 2 ' 

2a - 2«i + 2a 2 4a - 4«i + 4a 2 

-2ao + 2aj — 2a 2 — 2uq + 2aj — 2a 2 — 6ao + 3oti — 2a 2 
4«o — *i 4ao — 

— 4a + «i — a 



It is typographically impractical to depict solutions for n > 6 

To complete our solution to the original matrix equation (O, we follow [ 1 , p. 55] 
in noting that B satisfies © if and only if the matrix C = RBR, where R denotes 
the reversed identity matrix, satisfies the corresponding equation ([6). Thus the 
dimension of the solution space to our original equation |(2)isn + 2|~|l whenever 
p = (-!)». 
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